Vortex-antivortex dynamics and field polarity dependent flux creep in hybrid 
superconductor / ferromagnet nanostructures 
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Vortex-antivortex arrays (VAA) that are created in a type-II superconducting film by lattices 
of ferromagnetic dots with perpendicular anisotropy have been investigated. The highest critical 
current is shifted to a non-zero value of the applied field, and current - voltage characteristics show 
different regimes of vortex motion even in zero applied field due to the presence of the VAA. Creep 
of interstitial vortices is observed at low driving forces. This flux creep is strongly field polarity 
dependent. 

PACS numbers: 74.25.Qt 74.25. Sv 74.78. Fk 74.78.Na 



In the mixed state magnetic field penetrates type-II 
superconductors in the form of vortices carrying a quan- 
tized magnetic flux, which normally corresponds to a 
single flux quantum Q]. In order to avoid dissipation 
in superconductors it is necessary to prevent vortex mo- 
tion. Effective pinning can for instance be achieved by 
introducing columnar defects 13, 0, lattices of submi- 
cron holes (antidots) 0, 0, 0> or ferromagnetic dots 
[H El El- Besides enhanced critical currents, field 
polarity dependent flux pinning was observed by Mor- 
gan and Ketterson in superconducting Nb films with 
arrays of ferromagnetic Ni dots jlOJ , whose magneti- 
zation has an out-of-plane component. It was shown 
that in this system the critical current density j c de- 
pends on the mutual orientation of applied field and 
magnetization of the dots. Stable vortex-antivortex ar- 
rays (VAA) are predicted to exist in hybrid supercon- 
ductor/ferromagnetic systems even in zero external field 
due to the inhomogeneous stray fields of the ferromagnets 
[H H H E Sill El E M ED, . Despite strong 
interest of theory in this subject, experimental reports on 
VAA have been very rare |23j • Experiments on magnetic 
dots with out-of-plane magnetization have been carried 
out only for relatively weak stray fields in the sense that 
the stray field was screened by Meissner currents in the 
superconductor pi El Ei f . 

In this paper we will show that the vortex dynamics in 
type-II superconducting Pb films is strongly influenced 
by the presence of ferromagnetic Co/Pd dots with per- 
pendicular anisotropy. We will characterize the dynamic 
behavior of stray field-induced VAA by means of current 
(I) - voltage (V) characteristics and demonstrate pro- 
nounced field polarity dependence of the flux creep. A 
clear sign for the existence of VAA is the presence of dif- 
ferent regimes of vortex motion even in zero applied field. 
Fig. n shows a schematic drawing of the investigated 
sample. The structural and magnetic properties of this 
sample have been described elsewhere |27]]. The sam- 
ple consists of a type-II superconducting Pb film (thick- 
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FIG. 1: Schematic drawing of the investigated superconduct- 
ing / ferromagnetic hybrid structure. The superconductor is a 
85 nm Pb film, which is covered by an array of ferromagnetic 
Co/Pd dots with out-of-plane magnetization. 



ness tpb = 85 nm, critical temperature T c — 7.26 K) 
evaporated on a I nm Ge base layer on a Si substrate 
with amorphous Si02 top layer. The Pb film is cov- 
ered by a 10 nm Ge layer, which prevents the influence 
of the proximity effect with the metallic dot array. In 
order to avoid inhomogeneities of the current, this con- 
tinuous Gc/Pb/Ge trilayer is patterned into a transport 
bridge (width w = 200 /^m, distance between voltage 
contacts d = 630 /im) by optical lithography and chemi- 
cal wet etching. Measurements of the upper critical field 
of this transport bridge allow an estimation of the co- 
herence length £(0) = 34 nm and of the penetration 
depth A(0) = 49 nm. The transport bridge is covered 
by ferromagnetic dots using electron-beam evaporation 
and electron-beam lithography. The dots consist of a 
3.5 nm Pd base layer and a [Co(0.4 nm)/Pd(1.4 nm)1m 
multilayer with perpendicular magnetic anisotropy [28j. 
The dots are arranged in a square array with period 
L = 1.5 fiTa. They have a square shape with a side length 
of approximately 0.8 /im with slightly irregular edges. 
The vortex pinning properties of this hybrid sample were 
investigated by electrical transport measurements for two 
different magnetic domain states of the Co/Pd dots. Af- 
ter aligning the magnetic moments of the dots in positive 
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FIG. 2: (a) V(I) curves at T = 7.10 K in zero field after de- 
magnetization (m = 0) and after magnetization of the sample 
in a high positive field (m z > 0). (b) Field dependence of the 
critical current of the superconductor with the dots in the two 
different magnetic states. A voltage criterion of Vcr = 1 (iV 
was used. 



^-direction (ra z > 0) (by magnetizing the dots in a high 
positive applied field of hqH = +1 T perpendicular to the 
sample surface) the dots are in a single domain state |27j . 
After an ac demagnetization, the dots are in a multido- 
main state where the net magnetic moment of each dot m 
is zero (m = 0) 27]. Note that the fields applied during 
all measurement are much weaker than the coercive field 
of the dot array at room temperature \i§H coe = 150 mT. 
V(J) curves of the sample are measured in a Quan- 
tum Design Physical Properties Measurement System 
(PPMS) in standard four-point geometry by increasing I 
starting from 1 = 0. Figure Ufa) shows two V(I) curves 
measured with the sample in the m = and the m z > 
state at temperature T = 7.10 K in zero field. There is 
a remarkable difference between the two curves: In the 
magnetized state V shows a steep transition to the nor- 
mal state at about 6 mA, while in the m z > state this 
transition already appears around 1 mA. 
We define the critical current I c as the current where 



the voltage exceeds a certain voltage criterion V cr . The 
critical current density j c can be calculated using j c = 
Ic/{tpb w )- The field dependence of I c at T = 7.10 K 
is shown for V cr = 1 /iV in Fig. |2Jb) for both mag- 
netic states. The first matching field H\ of this sample 
corresponds to the field value which generates one flux 
quantum per unit cell of the dot array and is given by 
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bo/L 2 = 0.92 mT (</>o is the superconducting 



flux quantum). In the m = state the maximum of I c is 
observed at zero field, and the I C (H) curve is symmetric 
with respect to H = 0. Clear matching effects appear at 
positive and negative first matching field, indicating that 
the multidomain dots act as pinning sites for vortices due 
to the local suppression of the order parameter under the 
dots caused by the stray field, or due to the magnetic pin- 
ning of the vortices by the domain structure. 
In the m z > state I c is lower than in the m = state 
for H/H\ < 3.5. Surprisingly, the highest value of I c is 
not observed at H = 0, but at H/H\ = 2. Clear match- 
ing features are seen at H/Hi = 0, 1 and 2. The strong 
overall reduction of the critical current indicates that the 
shift is not due to a mere compensation effect of the ap- 
plied field by the local stray field of the dots |2jj , but that 
the vortex dynamics is strongly modified due to vortices 
induced by the dots, as will be explained below. 
The significant alteration of I C (H) which is obtained only 
by a different magnetic history of the sample must be re- 
lated to the changed magnetic domain state of the dots. 
It should be noted that H denotes only the externally 
applied field. The local stray field of the dots is also 
present at H = 0. However, when integrated over the 
superconducting layer, the net average stray field of the 
dots will be near zero because of the returning stray field 
of opposite sign. For the further analysis of the data we 
neglect the influence of heating effects, nonuniform cur- 
rent distributions and surface barriers since they are the 
same in both magnetic states and cannot account for the 
observed significant change in I C (H). In the following we 
will show that the existence of two vortices or one giant 
vortex with vorticity two directly beneath each dot can 
consistently explain the anomalous field-shift and also 
other features in the V(I) curves that will be discussed 
below. Additionally the returning stray field of the dots 
generates two antivortices per unit cell of the dot array 
in the superconductor. These antivortices are located at 
interstitial positions of the dot array. This assumption is 
justified by several theoretical calculations that predict 
this picture as a possible scenario [l3L H3. H5I l2l| . In the 
case of an isolated dot both, vortices and antivortices, 
are bound to the dot, but for a regular array of dots the 
periodically modulated stray field could change this pic- 
ture. If the spacing between the dots is sufficiently small, 
a motion of vortices and antivortices should be possible 
when a transport current exerts a sufficiently strong driv- 
ing force on them. 

Fig. shows the V(I) curves measured at H = 0, 
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FIG. 3: Field polarity dependent flux creep is observed in 
current - voltage characteristics of the superconductor at T = 
7.10 K. The dots are in the m z > state and V is plotted in 
logarithmic scale. The full lines are guides to the eye. The 
arrows indicate the crossover between low and high current 
regimes. The dotted lines indicate the used voltage criteria 
for determination of I C {H) (see text). The inset shows the 
field dependence of the activation energy Uo determined from 
the slope of the In V(I) curves. 

±0.7 mT, ±1.2 mT and ±1.7 mT with V plotted on 
logarithmic scale. Low and high current regimes must be 
distinguished: In the low current regime only interstitial 
flux lines move, whereas at higher applied current both 
interstitial flux lines and the ones that were originally 
pinned on the dot sites move. Different dynamic regimes 
for low and high currents have also been observed in V(I) 
curves of superconducting films with arrays of antidots 
|29) and arrays of magnetic dots (sjj. Thus, the V(I) 
curve at H — is in agreement with the assumption that 
interstitial antivortices are present even in zero field. 
In the low current region (V < 5 /iV) the InV(I) curves 
measured in zero field and positive fields have a linear 
slope as is indicated by the straight lines. This linear 
slope is less clear for the curve at —0.7 mT and disappears 
as H is further decreased. A linear \nV(I) dependence is 
predicted in the Anderson and Kim model for flux creep 
|31| . In this model the electric field E generated by the 
vortices creeping from one pinning site to another in the 
limit I rs I c is given by 

£oce^^^ I ~ Ic, (1) 

with ks the Boltzmann constant and Uq — Uq(T, H) the 
activation energy, i.e. the height of the potential barrier 
between two adjacent pinning sites. Thus, 

^V = ^-(y-l) + C{H,T), (2) 



where C is independent of I. This means that in zero 
applied field the experimentally observed ln(V) oc I be- 
havior is consistent with a creep of antivortices over the 
potential barriers that are formed by the dots and the 
induced vortices. Hence, Fig. [3] shows that flux creep 
in hybrid superconductor /ferromagnet systems is strongly 
depending on the field polarity. The exact creep pattern 
of the antivortices in the sample remains an open ques- 
tion. An antivortex can for example move by anihilating 
first with a vortex on the edge of a dot, and at the same 
time a new antivortex is generated at the opposite dot 
edge. Another possibility is that the antivortices creep 
through channels between the dots. A theoretical treat- 
ment in the framework of the time-dependent Ginzburg- 
Landau theory is required to answer this question. 
The slope of the In V(I) curves is only determined by Uo, 
T and I c . For the calculation of Uo (results are shown in 
the inset of Fig. O we have defined I c with a criterion of 
V cr = 0.2 jtiV. Up to H2, both Uo and I c increase with H , 
which is related to the decreasing number of antivortices. 
As the number of antivortices goes down, the effective 
barrier that they have to overcome is increased due to a 
decreasing mutual repulsive interaction. Note that this 
interaction is not included in the model (in its simplest 
interpretation). At Hi a peak in Uo{H) is observed, be- 
cause of the commensurability of the antivortices with 
the periodic pinning potential. Above H2, Uq starts to 
decrease drastically as H is increased. This crossover 
can be explained by the appearance of vortices at inter- 
stitial positions of the dot array: the number of vortices 
increases and therefore Uo decreases with increasing H . 
The most important conclusion of these results is that 
the observed features are in good agreement with the as- 
sumption that at H = two antivortices per unit cell are 
induced in the superconductor at interstitial positions of 
the dot array. Small inhomogeneities of the dots in a real 
sample could lead to a variation of the stray fields of the 
dots, which for a small fraction of dots could result in 
three or only one induced vortex-antivortex pair instead 
of two. However, this will not influence the general be- 
havior of the sample. 

A crossover to a linear V(I) dependence is observed for 
H/Hi > 2.8. The main panel in Fig. 0] shows the low 
voltage part of the V(I) curves in the m z > state for 
H/Hi = 2, 3, 4 and 5.2. The linear dependence is con- 
sistent with the situation that local pinning forces are 
overcome resulting in a steady motion of the interstitial 
vortices j3^| with velocity v given by 

r)v = (/ - I c )cj>o/{tp b w), (3) 

with r] a viscous drag coefficient, which can be determined 
from 

V = Ed = nf (j>ovd = n f (f>od(I - I c jiow)/(r) t Pb w), (4) 

where nt is the areal density of moving interstitial 
vortices given by nj — Ho\H — H2\/4>o- Ujiow is defined 
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FIG. 4: Current - voltage characteristics of the superconduc- 
tor at T = 7.10 K in the m z > state at HjH x = 2, 3, 4 
and 5.2 for small voltages. The slope of the straight lines de- 
termines the viscous drag coefficient rj in the flux flow regime 
H > 2.8. 

as the value of I where the straight lines characterizing 
the flux flow regime indicated in Fig. 01 cross the /-axis. 
From the slope of the V(I) curves and the value of 
Icfiow, 1S calculated using eq. |@J, see the inset of 
Fig. 0] In the Bardeen-Stephen model for free flux flow 
|33j it is assumed that r\ does not depend on the applied 
field. The flux flow phase of the interstitial vortices in 
our measurements can therefore not be interpreted in 
terms of free flux flow. The mutual interaction between 
flowing interstitial vortices and pinned vortices on the 
dot sites can not be neglected and should be the main 
reason for the increase of r) as H is increased. 
Finally, we would like to point out that ferromagnetic 
dots can be very useful for the new field of fluxon- 
ics. Here the aim is to guide the vortices through 
nanoengineered arrays of pinning sites or channels in 
the superconductor j^, which makes it possible to 
develop new devices, like superconducting reversible 
rectifiers |3(| or ratchet cellular automata based on 
nanostructured superconductors |34j . Using VAA in 
a well-designed fashion makes it possible to operate 
fluxonic devices without exposing them to externally 
applied fields. For instance, by using triangular shaped 
magnetic dots it is possible to measure DC voltages 
upon application of an AC current |3(ij ] even in zero field. 
Furthermore, by using rectangular dots that generate 
one vortex-antivortex pair per dot, logic elements can 
be developed [34| where the and 1 state are defined by 
the position of the antivortex. 

In conclusion, vortices in hybrid superconduc- 
tor/ferromagnet nanostructures can have peculiar 
dynamic properties due to the presence of the vortex- 



antivortex arrays that are induced in the superconductor 
by the inhomogeneous stray field of the ferromagnet. 
In zero applied field two regimes of vortex motion are 
observed, which proofs the existence of these induced 
vortex-antivortex arrays. For the first time field polarity 
dependent flux creep is demonstrated. 
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